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lO . Abstract 

O 

We study Uranium isotopes and surrounding elements at very large neu- 

r^ ■ tron number excess. Relativistic mean field and Skyrme-type approaches 

with different parametrizations are adopted in the study. Most models show 

clear indications for isotopes that are stable with respect to neutron emission 

far beyond N=184 up to the range of around N=258. 

/\ • Keywords: neutron-rich nuclei, drip line, uranium isotopes 

1. Introduction 

One of the foremost topics in modern nuclear physics is the study of 
exotic nuclei far beyond the stability line. A wide range of current experi- 
mental activities as well as new upcoming facilities, like FAIR at GSI and 
FRIB at MSU, will expand the knowledge of radioactive nuclei up to or close 
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to the neutron drip line. A key issue of these investigations is to deepen 
the understanding of the nuclear interactions at large isospin asymmetries, 
which could also yield important input for the study of neutron stars, as the 
most neutron-rich system known. The knowledge of the properties, espe- 
cially binding energies and estimates for lifetimes of neutron-rich isotopes is 
also crucial input for r-process nucleosynthesis calculations. In recent years 
there have been various studies of the structure of the neutron drip line 
Q, 3 3 i, B B 0, B B E [ll|. in the latter publication (ll| particle-stable 



Uranium isotopes up to a neutron number of 204 were calculated. A number 
of recent studies suggest a drip line that might be less smooth as originally 
expected. Results show potent ial p eninsulae of isotopes sticking out into the 



area of unstable nuclei J12l . Il3l Il4j . Il5| . Other calculations suggest islands of 



(meta-)stable nuclei beyond the drip line jl6|. 

In the following we study very heavy nuclei, from Uranium to Rutherfordium 
and investigate the final neutron number that can be bound by the nucleus. 
The results show nuclei with extremely large neutron excess that are still 
stable with respect to neutron emission. This behavior is seen in a number 
of model approaches and parameter sets. 

2. Theoretical Approach 

The calculational framework used is a relativistic mean-field model (RMF) 
with a parametrization fitted to binding energies, diffraction and rnis radii, 
and surface thicknesses of a set of nuclei jlTj as well as non-relativistic 
Hartree-Fock Skyrme-type models for various sets of Skyrme forces that were 
also used in the study of neutron-rich nuclei. 

The calculations were performed in spherical symmetry as well as in runs 
with deformation assuming axial symmetry. The equations were solved on a 
spatial grid in one and two dimensions to model the nuclei. In comparison 
the Skyrme calculations were also done in a self-adapting harmonic oscillator 
basis, where the oscillator asymmetries are adjusted iteratively to minimize 
the binding energy by allowing for nuclear deformations. 

We employed a BCS-type pairing, where we used delta force pairing for the 
relativistic and Skyrme calculations on the grid and constant force for the 
calculations using a harmonic oscillator basis. Certainly, improvements of 
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Figure 1: Total binding energy of Uranium isotopes for different parameterizations. The 
calculation was done for spherical nuclei. The curves show the result for a relativistic 
mean-field (NL-Z2) ^ and Skyrme models (SkM*, SkM, SLy6). as determined on a 
spatial grid. 



the pairing description should be done in future investigations, however, the 
basic quahtative results most likely will not be affected. 



3. Numerical Results 

Our numerical results of a spherical calculation of Uranium isotopes are 
shown in Fig. 1. One can see the total energy of the nuclei as function 
of neutron number. The curve exhibits an absolute minimum at ^^^JJ (in 
the case of the relativistic force), signaling the last neutron-stable isotope. 
Several results for RMF and Skyrme calculations as well as for a relativistic 
effective chiral Lagrangian jl8|] have been obtained most of them exhibiting 



Model N(ld) N(2d) Bad [MeV] 



NL-Z2 [17 


258 


258 


2021.34 


SkM* [19| 


246 


220 


1987.39 


SkI4[20] 


258 


218 


1975.41 


SLy6[2i| 


184 


206 


1928.92 


Xm118| 


184 


184 


1957.56 



Table 1: Heaviest Uranium isotopes for Id and 2d calculations. 

a mininiuin at large neutron number. In the cliiral Lagrangian the bary- 
onic masses are generated via spontaneous symmetry breaking that generate 
non-vanishing vacuum expectation values of the scalar fields. The quality 
of the model in reproducing knwon nuclear structure data is comparable to 
other relativistic models [18']. As it is based on a fiavor-SU(3) description 
the nucleons also couple weakly to the strange scalar field that corresponds 
to the strange quark condensate. This approach does not show a significant 
increase of bound neutrons with Z. Here a shell closure around A^ = 256 only 
develops at much higher atomic charge oi Z > 116. 

In order to check for possible effects from the grid size of the calculation 
we repeated the calculation for various box sizes from 30 fm to 80 fm in the 
case of Uranium without observing significant changes of the binding energy 
beyond the range of standard numerical uncertainties. 



Fig. [2] shows the result of a two-dimensional calculation including de- 
formation assuming axial symmetry. The calculations were performed on a 
spatial grid using three different initial deformations (oblate, spherical, and 
prolate) for each nucleus. The relativistic force maintains its most bound 
nucleus for large neutron number N=258, whereas the Skyrnie forces show 
shifted minima as can be observed in the figure. For comparison. Table [1] 
lists the heaviest uranium isotopes for the different parameter sets in spher- 
ical approximation and two-dimensional calculations. 

The observed behavior is not restricted to Uranium but hold true for sur- 
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Figure 2: Total binding energy of Uranium isotopes for different parameterizations. The 
calculation was performed for axially symmetric deformed nuclei on a two-dimensional 
spatial grid. The curves show the result for a relativistic mean-field (NL-Z2) and Skyrme 
models (SkM*, SkM, SLy6). 



NL-Z2 SkM* SkM SLy6 xm 



92 


258 


220 


218 


206 


184 


94 


258 


230 


230 


208 


184 


96 


258 


258 


230 


208 


202 


98 


258 


258 


258 


218 


208 


100 


258 


258 


258 


220 


210 


102 


258 


258 


258 


220 


214 


104 


258 


258 


258 


232 


220 



Table 2: Isotopes with maximum binding energy for different models. 
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Figure 3: Single-particle levels of neutrons for Uranium isotopes assuming spherical sym- 
metry using the NL-Z2 parametrization. The line shows the last occupied HF orbit ne- 
glecting pairing. The main cause of the semi-shell closure originates from the filling of the 
Ik ^+ state. The other bound states shown are (in the order as seen for N=220, starting 



from the lowest level) 2h 



3f 



4p §-, 4p 



3f 



bound for higher values of N as can be seen in the plot. 



The Ik ^- 



state is only 
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Figure 4: Deformation P2 for different Californium isotopes for the parameter set SkM'' 



rounding elements. Table [2] shows the most bound isotopes for a range of 
Z between 92 and 104. Those isotopes exhibit vanishing quadrupole moment. 

The level diagram of the neutron levels from 2.5 MeV below the contin- 
uum for the case of Uranium isotopes is shown in Fig. |3] . The lfcl7/2"'" 
level, which has a high degeneracy, drops with neutron number effectively 
creating a new shell for large total neutron number of A^ = 258. A similar 
effect occurs in the case of (some) Skyrme forces as seen in Table |2j 

As an example of the isotope behavior of the nuclear deformation. Fig. 
m shows the axial deformation as it changes towards the spherical state at 
A^ = 258 for the parameter set SkM* for Californium isotopes. The binding 
energy as function of the axial deformation for three isotopes along this iso- 
tope chain in a calculation including a deformation constraint is depicted in 
Fig. |5l Here one can observe the the change from prolate to oblate, and to 
finally a spherical state. 

Looking at the density distribution of protons and neutrons in the case of the 
extreme nucleus ^^^JJ (NL-Z2) one observes an outer shell of essentially pure 
neutron matter as can be seen in Fig. [61 The neutron skin of this nucleus 
results as 

Tskin = rn-rp = 7.72 fm - 6.64 fm = 1.08 fm. (1) 

The isospin asymmetry of such a nucleus 

is quite substantial compared to a value of 0.7 - 1 inside of a neutron star or 
a value of less than 0.22 for a ^°^Pb nucleus. 

Fig|7] shows the one and two-neutron separation energies Sn and S2n at the 
magic neutron number of A^ = 258. Here, S'„ was estimated using Koopmans' 
theorem |[22i]. The results for a spherical and a two-dinemsional calculation 
are shown, which are in agreement with each other pointing to the fact that 
at this neutron number the nuclei are spherical. One can observe that for 
Z = 92 S2n just crosses the zero line yielding nuclei that are stable with 
respect to 2 (and 1) neutron emission. However, as can be seen in Tabled 
for Z = 96 this state becomes the most stable nucleus in the isospin chain 
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Figure 5: Energy of different Cf isotopes with N = 222, 242, and 256 as function of axial 
deformation /32 . One can observe the change in deformation from prolate to slightly oblate 
to spherical. 
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Figure 6: Density distributions r^pB of protons and neutrons in ^so^j j-qj. ^j-^g parameter 
set NL-Z2. 
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Figure 7: One- neutron {Sn) and two- neutron [S^n) separation energies for neutron number 
N = 184 as function of Z for Skyrme parameters SkM*. The results were obtained in a 
spherical (SPH) and two-dimensional code (DEF) using a basis of oscillator states. 



for this particular Skyrme parametrization. Thus one can observe the devel- 
opment of similar metastable peninsulas along the neutron drip line as had 
been observed in 



4. Conclusions 

We have calculated heavy nuclei with extreme neutron numbers. Nuclei 
with up to 258 neutrons are shown to be stable in a range of self-consistent 
model calculations. It is interesting to note that one can achieve self-bound 
largely neutron system with an asymmetry of 0.48. It might be worthwhile 
to correlate the existence of the states in various models with neutron star 
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observables. At the moment the calculation is largely a theoretical exercise as 
it is not clear how to produce such nuclei in a realistic experimental setup in 
measurable quantities, even if, in principle, a lead-lead collision, for example, 
would contain enough neutrons. Recent speculations on using high-intensity 
pulsed neutron sources from reactors for producing large-N nuclei |2J| sound 
like a potentially very interesting new experimental idea for producing very 
neutron-rich nuclei, but this approach is still very far from practical real- 
ization. In the same vein, for now at least there is no scenario predicting 
production of such states in a realistic astrophysical environment, although 
similar states might occur in the crust of neutron stars. However, the knowl- 
edge of the boundary of particle-stable nuclei in itself is still a question of 
fundamental interest in studies of strong-interaction physics. 
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